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Correlating off-stoichiometric doping with nanoscale electronic disorder and 
quasiparticle interference pattern in high-Tc superconductor Bi2Sr2CaCu208+5 
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A microscopic theory is presented for the observed electronic disorder in superconducting 
Bi2Sr2CaCu208+i. The essential phenomenology is shown to be consistent with the existence of two 
types of interstitial oxygen dopants: those serving primarily as charge reserviors and those close to 
the apical plane contributing both carriers and electrostatic potential to the Cu02 plane. The non- 
linear screening of the latter produces nanoscale variations in the doped hole concentration, leading 
to electronic inhomogeneity. Based on an unrestricted Gutzwiller approximation of the extended t-J 
model, we provide a consistent explanation of the correlation between the observed dopant location 
and the pairing gap and its spatial evolutions. We show that the oxygen dopants are the primary 
cause of both the pairing gap disorder and the quasiparticle interference pattern. 

PACS numbers: 71.27.+a, 71.18.+y, 74.25.Jb, 74.70.-b 
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Remarkable electronic inhomogeneities have been ob- 
served by scanning tunneling microscopy (STM) in high- 
Tc superconductors BizSraCaCuzOs+^P, S Si, 0, 
Caa-^Na^CuOaClaBBl, and Bi2_:rPbxSr2CuO,j'8f over 
a wide range of doping. The hallmark of the inhomo- 
geneity is the disordered, nanometer scale variation of 
the pairing energy gap and its anti-correlation with that 
of the coherence peak height in the local density of states 
(LDOS). The origin of the electronic disorder has been 
the focus of several theoretical studies [9|, |lO|, |ll|, [ij, |l3| . 
It was emphasized [l|, |9|, [lOj that off-stoichiometric dop- 
ing a Mott insulator, such as the high-Tc cuprates, cre- 
ates inherently inhomogeneous electronic states associ- 
ated with the interstitial or substitutional dopant atoms. 
Ionized, these off-plane dopants act as disordered cen- 
ters of nonlinearly screened electrostatic potential, lead- 
ing to a spatially inhomogeneous distribution of the local 
doping concentration (LDC) and hence that of the local 
pairing gap. In a short coherence length superconductor 
described by the short-range resonance valence bond the- 
ory [l5|, lld l. the local pairing gap is anti-correlated with 
theLDC~P,[i3]. 

The notion of dopant induced electronic disorder is 
supported and further advanced by the recent STM ex- 
periments of McElroy et.al. [J|. Identifying the spectral 
peak around — 0.96eV in the LDOS with the presence of 
a local excess oxygen atom, the correlation between the 
dopant location and the gap inhomogeneity has been es- 
tablished ubiquitously. However, the observed dopants 
reside close to the regionos of large pairing gap. This 
is counter-intuitive since the negatively charged oxygen 
ions are expected to attract nearby holes and create a 
higher LDC with a smaller pairing gap [Ol, llO|. This dis- 
crepancy promotes the idea that dopant induced local 
structural distortions ma y p lay a more important role 
than potential disorder [J, ll7| and phenomenological the- 
ories in which the dopants serve as large pairing centers 
[l3| . The origin of the dopant induced electronic disorder 



has remained a central unresolved issue. 

In this paper, we show that, while dopant induced 
structural effects are present, the main cause for the 
electronic disorder is the dopant potential induced LDC 
modulations. A clue as to how the latter can be con- 
sistent with the dopant locations comes from the follow- 
ing observations made on the experiments of McElroy 
et.al. [J|: (1) The observed dopants associated with the 
— 0.96eV peak cannot account for the total number of 
oxygen dopants. For example, in the underdoped sample 
with A = 65meV, the observed dopant density is 2.7%. 
Even if every oxygen dopant is fully ionized, the average 
doping would only be 5.5%, which is much smaller than 
the expected value around 11% p^]. Thus, a substantial 
number of the dopants has not been located, too large to 
be accounted for by the uncertainty in the doping pro- 
cess. (2) The correlation between the observed dopants 
and the pairing gap is weak and well defined only in the 
statistical sense. An appreciable number of them can be 
found in the smaller gap regions or to straddle the bound- 
aries between the small and large gap regions. Thus, it is 
unlikely that the observed dopants strongly and directly 
affect the local electronic structure. (3) The in-gap, low 
energy states in the LDOS have their weight concentrated 
in regions away from the identified dopants. This sug- 
gests that these electronic states which are encoded with 
the quasiparticle interference modulations [J] are likely 
localized or pinned by additional confining potentials in 
regions away from the identified dopants. 

Based on these observations, we conjecture that there 
are two types (A and B) of interstitial oxygen dopants 
in Bi2Sr2CaCu208+5. The type-B dopants serve pri- 
marily as charge reserviors. They only couple weakly 
to the Cu02 plane. These near nonbonding oxygen or- 
bitals give rise to the — 0.96eV peak in the LDOS \^. 
Recent ARPES measurements show that the dopant in- 
duced states near — 0.96eV have Bi symmetry and do not 
mix with the doped holes residing in the planar orbitals 



|18l | . This further supports the identification of the for- 
mer with the type B-dopants. The type-A dopants, on 
the other hand, strongly affect the local electronic struc- 
ture in the Cu02 plane. Their electrostatic potential en- 
hances the LDC which in turn pushes the orbital energy 
of type-A dopants into the broad valence band spectra 
below — 1.2eV not yet accessible by STM probes. One 
possibility is that the type-B dopants sit above the BiO 
plane while type-A dopants come close to the apical SrO 
layer. We will discuss the microscopic origin of the latter 
at the end. 

To support this physical picture, we extend the t-J 
model to include the oxygen dopant potential, 

H = - Y^ UjPcl^Cj^P + JY^ (Si • Sj - -n^nj) 

i^j i 

Here c|^ creates an electron that hops between near 
neighbors of the Cu square lattice via tij . Repeated spin 
indices are summed and Hi = c\^Cin is the density opera- 
tor. The LDC is given by Xi = 1 — rii, Ui = {cl^Cic}- The 
average doping will be denoted as S = (l/Ng) J2i ^i on a 
lattice of Ng sites. The second line in Eq. ([T|) describes 
the long-range Coulomb interactions between the elec- 
trons in the plane and between the in-plane doped holes 
and the two- types of off- plane dopants, Vfj — Vc/\ri—rj\ 
and 
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where (Ia, ds are the setback distances and Na, Nb the 
number of type-A and type-B dopants respectively, 5 — 
2{Na + Nb)/N,. 

Eq. ^ describes doped Mott insulators because of the 
projection operator P that removes double occupation. 
The projection is most conveniently implemented using 
the Gutzwiller approximation by the statistical weighting 
factors multiplying the coherent states, thus renormaliz- 
ing the hopping and the exchange parameters [20| . Since 
the dopants break translation symmetry, we extend the 
approach to the disordered case by the renormalization 
tij -^ gljtij and J -^ g^J where the Gutzwiller factors 
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(3) 



depend on the local doping at the sites connected by 
the hopping and the exchange processes. The exchange 
term is decoupled in terms of the bond xtj = {na'^j'y) 



and the pairing A^ 



{^aa'CiaCja') ficlds, leading to a 



renormalized mean-field Hamiltonian, 



Hga= - ^9ljUjc\„Cj„ + 



i=^3 



7 , SjCio-Cio- 



~ 4^^ E Sij {I^Xj^aa'CiaCja' + XycLcjT + h.c] 

+ ij^4(|A,,f + |x.,f)-5]A.n,;. (4) 

The local energy for the electrons is given by Ei = Vsc (i) + 
\i — M/j where /x/ is the chemical potential and \i is a 
fugacity that, together with the last term, ensures the 
equilibrium condition under local occupation 
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21|. V^c(*) is the screened Coulomb potential. 






Y,,[i) = Y^^ + YP^Y, 



(5) 



This is the driving force of the electronic disorder through 
local doping variations. The correlation of the latter to 
the local pairing gap disorder is caused by the Gutzwiller 
factor g* in Eq. ([3]) that modulates the kinetic energy lo- 
cally. Minimizing the ground state energy of Eq. ([4]), we 
obtain the self-consistent equations for the set of param- 
eters (Ajj, Xij, Xi,Xi, Hf), which are solved iteratively for 
a given average doping S. 

We present our results for systems of 32 x 32 sites. 
We use tij = (0.48, -0.16, 0.05, 0.05, -0.05)eV and J = 
O.OSeV such that the quasiparticle dispersion in Eq. ((4]) 
agrees with that measured by angle-resolved photoemis- 
sion [l9|. For simplicity, half of the dopants are taken 
to be B-type and the other half A-type distributed ran- 
domly with dA — la. The bare electrostatic poten- 
tials are set to Ya ^ Vc = 0.5eV, and Yb = for 
the weak coupling between B-dopants and doped holes. 
Note that the locations of type-A and type-B dopants 
are naturally anticorrelated over the average interdopant 
distance a^/2jS. The spatial distribution of the LDC 
Xi and the d-wave pairing order parameter Ad{i) are 
shown as 2D maps in Figs, la and lb for a typical sys- 
tem at S = 10.2%. The projected dopant locations are 
superimposed. The A-dopants are ubiquitously corre- 
lated with the LDC. Due to the local pairing nature of 
the short-range RVB state, Xi is in turn strongly anti- 
correlated with Ad{i) [9|. Notice that the modulations 
induced by the A-dopants leaves the B-dopants unwit- 
tingly correlated overall with the low doping and strong 
pairing regions. The LDOS is calculated from the pro- 
jected retarded Green's function in the Gutzwiller ap- 
proximation 21], LDOS(i,w) = -2lmgliG{i,i,Lu + iO+). 
To reduce finite size effects, we average over different 
boundary conditions corresponding to 20 x 20 supercells. 
As in the STM experiments, the local tunneling gap At 
is extracted from the coherence peak position at positive 
bias in the LDOS. In Fig. Ic, the tunneling gap map is 
shown with the dopant locations. Evidently, At is small 
near the A-dopants where the LDC is high. Notice how- 
ever, that the B-dopants are found with high statistics in 
and around the regions where the tunneling gap is large, 
consistent with their identification with those observed 
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FIG. 1: Doping and pairing disorder on a 32 x 32 system 
at 5 = 10.2%. 2D maps are shown for the LDC Xi (a), di- 
mensionless d-wave pairing order parameter A^j (b), and tun- 
neling gap At in meV (c) with the projected type-A (open 
black circle) and type-B (white solid circle) dopant locations 
superimposed, (d) Correlation functions among the dopant 
location, At, Xi, and A^. 



by STM JJ. To elucidate the correlation between the 
dopants and electronic disorder, we calculate the nor- 
malized cross-correlation function between the dopants 
and the tunneling gap. 
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j:ds^T{nWEA^OAiB){r,: 



where Sf{ri) — f{ri) — (/). The dopant locations are 
modeled by Lorentzians of width 2a 4|. Fig. Id shows 
that, while At is strongly anticorrelated with A-dopant 
locations (Cat-Oa(O) ~ —0-6), the B-dopants (observed 
by STM) are positively correlated with the local tun- 
neling gap, and moreover, the correlation is significantly 
weaker, C'at=Pa(0) ~ 0.3, in excellent agreement with 
experiments y. 

The large LDC variation is not inconsistent with the 
STM integrated LDOS variation that diminishes only 
when integrated up to — 0.6eV to — 0.9eV, where large 
incoherent background and the spectral weight associ- 
ated with the dopants contribute appreciably [J|. It is a 
subtle task to relate quantitatively the integrated LDOS 
or the topography to the LDC since the former are ob- 
tained in the constant tunneling current mode where the 
tip to sample distance changes significantly, leading to 
reduced spatial variations of the spectral weight |l|, |9| . 

This form of dopant induced electronic disorder can de- 
scribe the basic properties of the inhomogeneous low en- 
ergy states observed experimentally. In Fig. 2, we present 
the LDOS along the two line-cuts indicated in Fig. Ic. 
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FIG. 2: Evolutions of the LDOS along two line cuts (shown 
in Fig. Ic) from small to average (a) and average to large (b) 
gap regions. Dashedlines are guide to eye. 



Clearly, larger (smaller) gap regions are associated with 
smaller (larger) coherence peaks in agreement with STM 
[l|, y, y, y, Isj . Note that the line-cut passes through lo- 
cally highly overdoped regions with Xi > 25%. In a uni- 
form system, the pairing gap would be vanishingly small 
and the system essentially in the normal state at such 
high doping levels. Fig. 2a and Fig. Ic show, however, 
that the smaller local gap is still sizable. That overdoped 
regions of sizes not exceeding the coherence length have 
gap values considerably above those in a clean sample at 
the corresponding dopings was first pointed out in Ref. [9| , 
and is likely a manifestation of the proximity effect due to 
the surrounding larger gap regions. Indeed, Fig. 2 shows 
that as we move away from the small (large) gap region, 
a larger (smaller) gap emerges and evolves progressively 
stronger until it becomes the dominant gap as the line 
cut enters the lower (higher) doping region. This feature 
has been observed recently by STM with high energy res- 
olutions y. 

Next we show the A-dopants, the culprit of strong gap 
disorder, also serve as the primary cause of the low- 
energy quasiparticle interference modulations [22]. In 
Fig. 3a, the 2D maps of the LDOS(z, us) is plotted at fixed 
energies to, showing different interference patterns. The 
maxima of the spatial modulations are preferentially cen- 
tered around the A-dopant sites, thus away from the B- 
dopants, consistent with STM observations. This turns 
out to be true for all low energy states. In Fig. 3c, we plot 
the local correlation between the LDOS and the dopant 
locations as a function of bias voltage. For both positive 
and negative energies, the correlations are positive and 
strong (reaching 0.9) with A-dopants and negative (an- 
ticorrelated) and weak (reaching —0.3) with B-dopants. 
The peak and dip around ±70meV are related to the Van 
Hove- like feature in the LDOS spectra seen in Fig. 2. Just 
as in the case of dopant-gap correlations, the strong cor- 
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FIG. 3: Electron LDOS maps (a) and the Fourier transform 
of the quasiparticle interference pattern (b) at low energies, 
(c) Local correlations of LDOS (i) and the dopant locations as 
a function of bias energy. 



relation of the interference modulations with A-dopants 
produces a weak anticorrelation with the B-dopants as in 
STM experiments [J]. To further illustrate the interfer- 
ence pattern, the Fourier transform of the quasiparticle 
LDOS are shown in Fig. 3b at the corresponding ener- 
gies. The dominant interference wavevector qi (marked 
by arrows) connecting the tips of the Fermi arcs [23| are 
clearly seen to disperse with energy near (0, ±27r/4a) and 
(±27r/4a, 0), while the qy connecting two tips of the same 
arc shows the opposite dispersion with e nergy , in remark- 
able agreement with STM observations 14l.l22l|. 

We now provide a possible microscopic origin of the 
type-A dopants. It is well known that Bi-based cuprates 
have a natural tendency toward Bi:Sr nonstoichiometry, 
i.e. a fraction of the Bi comes to the SrO apical plane 
and replaces the Sr in order to form the crystal structure 
[2i | . This creates the so-called A-site disorder [13, Ha] ■ 
Since the trivalent Bi'^^ replacing Sr^+ creates an ex- 
cess positive charge locally, it naturally attracts the neg- 
atively charged interstitial oxygen dopants to the vicinity 
of the apical plane, forming the A-dopants. The doping 
levels derived from the observed type-B dopants alone 
are about 5 ~ 6% less than the expected values in the 
samples studied by McElroy, et.al. [J|, suggesting about 
5 ~ 6% of the doped holes must come from the ^ 3% 
"missing" A-dopants. This is reasonably consistent with 
the typical Bi:Sr nonstoichiometry in Bi2212 where about 
5% of the Sr in each apical plane is replaced by Bi 24 1 . 
Recently, the effect of A-site disorder on Tc has been 



studied systematically by controlled trivalent (Ln^+) sub- 
stitution of Sr^+ in Bi2Sr2-i,Ln„CuOfi-i-A where Ln=La, 
Pr, Nd, Sm, Eu, Gd, and Bi [13, [ll]. We propose that 
the STM experiments be carried out on Bi2212 samples 
with controlled trivalent substitution of Sr in the apical 
plane. Our theory predicts that the density of the ob- 
servable type-B dopants would decrease with increasing 
Ln'^"'" substitution, while that of the type-A dopants in- 
creases, leading to stronger electronic disorder. 

We have shown that the electrostatic potential of the 
off-stoichiometric A-dopants can be the primary cause 
of the electronic disorder and quasiparticle interference 
modulations observed in Bi2Sr2CaCu208-i-5. The elec- 
tronic inhomogeneity in our theory is driven by that of 
the kinetic energy or the coherence of doped holes in a 
doped Mott insulator. Incoherent excitations beyond the 
Gutzwiller approximation and the dopant induced struc- 
tural distortions can also contribute to the electronic dis- 
order. We expect such dopant induced electronic disorder 
in all doped cuprate superconductors, with specific prop- 
erties dependent on the dopant locations and the crystal 
field environment, interstitial or substitutional, ordered 
or disordered. The off-plane dopant structure together 
with the role of the apical oxygen may account for the 
varying properties of the cuprates that share, otherwise, 
identical Cu02 planes. 
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